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ABSTRACT 

We present discovery imaging and spectroscopy for nine new z ~ 6 quasars found in the Canada-France High- 
z Quasar Survey (CFHQS) bringing the total number of CFHQS quasars to 19. By combining the CFHQS with 
the more luminous SDSS sample we are able to derive the quasar luminosity function from a sample of 40 
quasars at redshifts 5.74 < z < 6.42. Our binned luminosity function shows a slightly lower normalisation and 
flatter slope than found in previous work. The binned data also suggest a break in the luminosity function at 
Mi45() « -25. A double power law maximum likelihood fit to the data is consistent with the binned results. The 
luminosity function is strongly constrained (1 a uncertainty < 0.1 dex) over the range -27.5 < M1450 < -24.7. 
The best-fit parameters are ^(Mi^^q) = 1.14 x lO'^^Mpc'-'mag"', break magnitude Mi^^q = -25.13 and bright 
end slope /3 = -2.81. However the covariance between /3 and MJ'^jy prevents strong constraints being placed 
on either parameter. For a break magnitude in the range -26 < M^^^q < -24 we find -3.8 < /? < -2.3 at 95% 
confidence. We calculate the z = 6 quasar intergalactic ionizing flux and show it is between 20 and 100 times 
lower than that necessary for reionization. Finally, we use the luminosity function to predict how many higher 
redshift quasars may be discovered in future near-IR imaging surveys. 
Subject headings: cosmology: observations — quasars: general — quasars: emission Unes 



1. INTRODUCTION 

Observations of the most distant quasars provide important 
information on the state of the universe within the first billion 
years. They are used to probe the transition from a mostly 
neutral to mostly ionized intergalactic medium (IGM). They 
provide the only direct estimates of the metallicities and evo- 
lutionary states of massive galaxies. Quasars also probe the 
very early growth of supermassive black holes in the centers 
of massive galaxies. 

The success of surveys such as the Sloan Digital Sky Sur- 
vey (SDSS; Fan et al. 2006a; Jiang et al. 2009) and Canada- 
France High-z Quasar Survey (CFHQS; Willott et al. 2009) 
means that the number of known quasars at z ^ 6 is approach- 
ing 50. These samples are now large enough to determine 
the typical properties of quasars over a range of luminosities. 
One of the most remarkable results is that in most respects, 
quasars at z 6 have properties such as metallicity, emission 
line strength, hot dust mid-IR luminosity and cool dust far-IR 
luminosity very similar to lower redshift quasars (Freudling 
et al. 2003; Iwamuro et al. 2004; Fan et al. 2004; Jiang et al. 
2006; Wang et al. 2008). 

The number of known z 6 quasars is now large enough to 
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allow the luminosity function to be determined. The luminos- 
ity function is important since it encodes information about 
the build-up of supermassive black holes and gives the hard 
ionizing radiation output of the quasar population. Early work 
on the z 6 quasar luminosity function was based on deter- 
mining the normalisation, $, and bright-end power law slope, 
P, from the SDSS. This was done by direct determination 
from the known SDSS quasars by Fan et al. (2004) yielding $ 
a factor of 30 times lower than at z = 3 and /3 = -3.2 (95% con- 
fidence range of -4.2 to -2.2). Further constraints came from 
the lack of gravitational lenses amongst the known SDSS 
quasars which implied -4 (Fan et al. 2003; Richards 
et al. 2004). Jiang et al. (2008; 2009) combined the SDSS 
main sample with the SDSS deep stripe to find a somewhat 
shallower slope of /3 = -2.6 ± 0.3. However, the binned data 
shown in Jiang et al. (2009) suggest that there may be a break 
in the luminosity function at Mi4go « -26, in which case /? 
could be closer to the value determined by Fan et al. (2004). 

There have been several small area surveys searching for 
faint z 6 quasars. These are summarized by Shankar & 
Mathur (2007). Only the key results for surveys at z > 5.7 
are repeated here. Cool et al. (2006) discovered three quasars 
at z > 5, including one quasar at z = 5.85 with z' = 20.68, in 
the AGES survey which covers only 8.5 sq deg. As shown 
by Jiang et al. (2008), Cool et al. were very fortuitous to 
discover such a bright quasar at z > 5.7 in such a small sky 
area. The space density of similar luminosity quasars deter- 
mined by Jiang et al. from the much larger SDSS deep stripe 
is about six times lower than that of Cool et al. 

Mahabal et al. (2005) searched for high-redshift quasars 
around known SDSS quasars and serendipitously discovered 
one with magnitude z' = 23.0 at z = 5.70 foreground to a SDSS 
quasar located at z = 6.42. In total, Mahabal et al. searched 
2.5 sq deg. to a depth of z' = 24.5. Willott et al. (2005a) 
searched the first data release (TOOOl) of the CFHT Legacy 
Survey Deep and found no quasars to z' = 23.35 in 3.8 sq. 
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deg. The Mahabal et al. and Willott et al. results reveal that 
there must be a break in the luminosity function between these 
low luminosities (M1450 w -23) and the high luminosity SDSS 
quasars (M1450 « -27). Because low luminosity quasars pro- 
vide the bulk of ionizing photons from AGN, Willott et al. 
(2005a) showed that quasars were incapable of reionizing the 
IGM at z = 6, consistent with constraints from the unresolved 
X-ray background (Djikstra et al. 2004). 

In this paper we present the discovery of nine further 
CFHQS quasars and use a combined CFHQS/SDSS sample 
to derive the z = 6 quasar luminosity function. Sec. 2 presents 
the data on these quasars and notes on each one are given 
in Sec. 3. Sec. 4 discusses the depth, area and completeness 
of the parts of the CFHQS and SDSS which are used in this 
paper to determine the luminosity function. Sec. 5 presents 
binned and parametric derivations of the z = 6 quasar luminos- 
ity function. Sec. 6 discusses what this new derivation implies 
for the quasar intergalactic ionizing flux and search strategies 
for even higher redshift quasars. We present the conclusions 
in Sec. 7. 

All optical and near-IR magnitudes in this paper are 
on the AB system. Cosmological parameters of Hq = 
70 km s"' Mpc"', flm = 0.28 and = 0.72 (Komatsu et al. 
2009) are assumed throughout. 

2. NEW CFHQS QUASARS 

The optical imaging phase of the CFHQS and associated 
Canada-France Brown Dwarf Survey (CFBDS) is now com- 
plete. It draws on the Very Wide (VW), Wide and Deep 
components of the CFHT Legacy Survey (CFHTLS), the 
Subai-u/XMM-Newton Deep Survey (SXDS) and the Red- 
sequence Cluster Survey (RCS-2). For some of these surveys, 
additional imaging was obtained in PI mode to complete the 
necessary filter coverage. Further details of the observations 
are given in Willott et al. (2009). Near-IR and spectroscopic 
follow-up is still ongoing, so this paper does not present the fi- 
nal results of the CFHQS. Sec.|4]discusses the complete parts 
of the CFHQS which are used to determine the quasar lumi- 
nosity function in this paper. 

Quasar candidates were selected using the same color crite- 
ria as in Willott et al. (2009). The first cut is at i'-z' > 2 
to select very red objects which are a mixture of brown 
dwarfs and high-redshift quasars. Pointed /-band photometry 
with the ESO New Technology Telescope and Nordic Opti- 
cal Telescope is then used to separate these two types of ob- 
ject, the dwarfs having redder z' -J colors than the quasars. 
Spectroscopy of candidates was carried out mostly using the 
GMOS spectrographs at the twin Gemini telescopes. Two 
quasars (CFHQS J0050H-3445 and CFHQS J2229H-1457) were 
observed with the ESI spectrograph at Keck-II. Spectroscopic 
reductions were performed as detailed in Willott et al. (2007). 

Table[T]gives the positions and photometry for the nine new 
quasars. Table|2] details the spectroscopic observations. Note 
in this table and the rest of the paper, quasar names use an 
abbreviated form. Fig.[T]shows cutout images centred on each 
quasar and Fig.|2]plots their optical spectra. Larger area find- 
ing charts are in the Appendix. 

3. NOTES ON INDIVIDUAL QUASARS 

3.1. CFHQS J0050+3445 

The redshift of the quasar is z = 6.25 based on the narrow 
peak of Ly a emission on top of the broad line. This is consis- 
tent with the redshift based on the broad Mg II emission line 
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CFHQS J1429+5447 




CFHQS J2242+0334 

Figure 1. Images in the z' and J filters centered on the nine CFHQS 
quasars. Each image covers 20" X 20". The images are oriented with north 
up and east to the left. Two quasars are located close to the edges of the chips 
at i' and z' and therefore do not have data across the whole region. 
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Figure 2. Optical spectra of the nine newly discovered quasars. The expected locations of Ly a, Ly /3 and Ly 7 are marlced with dashed lines. All spectra are 
binned in 10 A pixels. The dashed Ly a lines often do not appear to line up with the peak flux in these plots (even though most redshifts are measured from the 
Ly a peak), due to a combination of the biiming and asymmetric line profiles. 
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Table 1 

Quasar positions and photometry 



Quasar 


RA and DEC (J2000.0) 


i' mag 


Z'l 


mag 


/mag 




-z' 


z'-J 


CFHQS J005006+344522 


00:50:06.67 


+34:45:22.6 


23.49 ±0.04 


20.47 


±0.03 


19.89 


±0.04 


3.02 


±0.05 


0.58 ±0.05 


CFHQSJOl 3603+022605 


01:36:03.17 


+02:26:05.7 


> 24.72^ 


22.10 


±0.09 


22.09 


±0.22 


> 


2.62 


0.01 ±0.24 


CFHQS J142952+544717 


14:29:52.17 


+54:47:17.7 


23.88 ±0.06 


21.45 


±0.03 


20.64 


±0.07 


2.43 


±0.07 


0.81 ±0.08 


CFHQS J022 122-080251 


02:21:22.71 


-08:02:51.5 


25.60 ±0.21 


22.63 


±0.05 


22.03 


±0.14 


2.97 


±0.22 


0.60 ±0.15 


CFHQS J222901+145709 


22:29:01.65 


+14:57:09.0 


24.89 ±0.16 


22.03 


±0.05 


21.95 


±0.07 


2.86 


±0.17 


0.08 ± 0.09 


CFHQS J210054-171522 


21:00:54.62 


-17:15:22.5 


24.55 ±0.21 


22.35 


±0.09 


21.42 


±0.10 


2.20 


±0.23 


0.93 ±0.13 


CFHQS J031649-134032 


03:16:49.87 


-13:40:32.3 


24.57 ±0.19 


21.72 


±0.08 


21.77 


±0.17 


2.85 


±0.21 


-0.05 ±0.19 


CFHQS J105928-090620 


10:59:28.61 


-09:06:20.4 


23.21 ±0.04 


20.82 


±0.03 


20.79 


±0.07 


2.39 


±0.05 


0.03 ±0.08 


CFHQS J224237+033421 


22:42:37.55 


+03:34:21.6 


24.22 ±0.12 


21.93 


±0.04 


22.13 


±0.12 


2.29 


±0.13 


-0.20 ±0.13 



Note. — All magnitudes are on the AB system. 
^ Where not detected at > 2a significance, a 2a lower limit is given. 



Table 2 

Optical spectroscopy observations of the new CFHQS quasars 



Quasar 


Redshift 


Date 


Resolving 


Slit Width 


Exp. Time 


Seeing 


Mi450 




z 




Power 


(Arcsec) 


(s) 


(Arcsec) 




CFHQS J0050+3445 


6.25 


2008 Sep 25 


4800 


1.0 


6250 


1.3 


-26.62 


CFHQSJOl 36+0226 


6.21 


2009 Aug 02 + 2009 Sep 13 


1300 


1.0 


5400 


0.9 


-24.40 


CFHQS J 1429+5447 


6.21 


2009 Mar 31 + 2009 Apr 30 


1300 


1.0 


5400 


0.9 


-25.85 


CFHQS J022 1-0802 


6.16 


2008 Oct 24 


1300 


1.0 


5400 


0.7 


-24.45 


CFHQS J2229+1457 


6.15 


2008 Sep 24 


4800 


1.0 


3600 


0.9 


-24.52 


CFHQS J2100-1715 


6.09 


2008 Jul 30 


1300 


1.0 


3600 


0.7 


-25.03 


CFHQS J0316-1340 


5.99 


2009 Aug 02 


1300 


1.0 


3600 


0.7 


-24.63 


CFHQS J1059-0906 


5.92 


2009 Mar 22 


1300 


1.0 


3600 


0.6 


-25.58 


CFHQS J2242+0334 


5.88 


2008 Oct 18 


1300 


1.0 


3600 


0.5 


-24.22 



Note. — GMOS spectra have resolving power 1300 and ESI spectra have resolving power 4800. As in Willott et al. (2009), absolute magnitudes (Mi 450) 
are calculated using the measured 7-band magnitudes and assuming a template quasar spectrum. Note that this is slightly different to the method used in 
Willott et al. (2007), which was based on measuring the continuum redward of Ly a in the spectrum and assuming a spectral index. 
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(Willott et al. in prep.)- This is the second highest redshift 
in the CFHQS. It is also the second most luminous CFHQS 
quasar with M^so = -26.62, which is nearly as luminous as 
some of the SDSS quasars of Fan et al. (2006b). 

This quasar is included in the group of RCS-2/CFHTLS 
VW quasars, however it is the only one in that group which 
is not actually located in these survey areas. We carried out 
z' band imaging of regions already observed at /' band with 
CFHT of the intragroup environment of the Local Group (PI 
R. Ibata). These observations were necessary to use sched- 
uled time (supposed to be /' band followup of RCS-2 fields) 
when the /' filter broke and was unusable. This field proved 
not to be ideal for the CFHQS because of several Local Group 
variable red giant stars which had non-contemporaneous i'z'J 
photometry appearing to be z ^ 6 quasars, which required re- 
peat photometry to eliminate. 

3.2. CFHQS JOl 36+0226 

The spectrum shows a fairly narrow broad Ly a line. The 
redshift of z = 6.21 is derived from the peak of the narrow 
Lya. This is the only quasar in this paper which was not 
detected in the /' filter. It has a lower limit on the color of /' - 
z' > 2.62. The z' -J = 0.01 color is quite blue for such a high 
redshift quasar This is likely because of the high equivalent 
width Ly a emission. 

The faint galaxy 3 arcsec north of the quasar (Fig.[T]i hap- 
pened to be located along the spectrograph slit. It has a single 
emission line at 7468 A which is likely [O ll] at z = 1 .00. The 
low luminosity of this galaxy (at observed frame wavelengths 
from /' to J band) indicates it has only a moderate mass and at 
this separation from the quasar will provide negligible gravi- 
tational lensing magnification (see e.g. Sec. 4 of Willott et al. 
2005b). 

3.3. CFHQS J1429+5447 

This quasar was discovered in the CFHTLS Wide W3 re- 
gion. However, it is much brighter than our Wide magnitude 
limit of z' = 23 and has the second brightest absolute magni- 
tude in this paper The spectrum shows a strong continuum 
with only a weak Ly a emission line at z = 6.21. This is con- 
sistent with the relatively red color of z' - / = 0.8 1 . The spec- 
trum shows a dark Lya absorption trough from 5.85 < z < 
6.16 with lots of sharp Ly a peaks at lower redshift. A more 
sensitive and higher spectral resolution spectrum is required 
for a complete analysis. 

There is an obvious Mgll absorption doublet at 9061.5, 
9083.8 A which corresponds to z = 2.241 . The apparent broad 
absorption at 9300 A is telluric absorption which has not been 
corrected for (and can be seen for several other quasars in 
Fig.©. 

3.4. CFHQS J022 1-0802 

Another quasar from the CFHTLS Wide (Wl field) and the 
faintest at z' band in this paper The Ly a emission peaks at 
z = 6.13 with a narrow spike atop an asymmetric broad com- 
ponent. There is also a likely N V A 1240 line consistent with 
the Ly a redshift. On the other hand, the near-IR spectrum 
(Willott et al. in prep.) shows the Mg II line at z = 6.16 which 
we adopt as the systemic redshift for this quasar 

3.5. CFHQS J2229+1457 

This is another CFHQS quasar with a strong, very narrow 
Ly a emission line. The best-fit redshift to the peak of the 



Ly a line gives z = 6.15. An unpublished near-IR spectrum 
shows the Mg II line to have an identical redshift. 

3.6. CFHQS J2100-17 15 

This quasar has a very similar spectrum to CFHQS J0221- 
0802 with a narrow Ly a peak at z = 6.09 (consistent with the 
unpublished Mg II line redshift). There is an interesting peak 
in the Lya forest at z = 6.00 which appears too far from the 
systemic redshift to be within the quasar-ionized near-zone 
for a quasar of such moderate luminosity. CFHQS J2 100- 
1715 has a rather red color of z' -J = 0.93 for its redshift, 
which may be partly due to the weak Ly a line, but also indi- 
cates a red spectral slope. 

3.7. CFHQS J03 16-1 340 

The spectrum shows a double-peaked Ly a emission line. 
We take the redshift of z = 5.99 from the blueward peak. In 
terms of the spectrum and photometry, there is nothing un- 
usual about this quasar. 

3.8. CFHQS Jl 059-0906 

This is one of the few CFHQS quasars to have a very broad 
Ly a line, typical of those in the SDSS main sample. The 
quasar is fairly bright (z' = 20.82), consistent with the idea 
that line width correlates with luminosity due to the black 
hole mass. There is no narrow peak to Ly a and consequently 
a large uncertainty on the redshift which is estimated to be 
z = 5.92 based on the full broad Lya profile and expected 
asymmetry due to hydrogen absorption. 

3.9. CFHQS J2242+0334 

Yet another quasar with a narrow Ly a peak just redward 
of a sharp cutoff. At z = 5.88 and M1450 = -24.22 this is the 
lowest redshift and least luminous quasar from the RCS-2 and 
CFHTLS Very Wide. It is also the bluest quasar in the CFHQS 
with z'- / = -0.2. This is not too surprising because the ex- 
pected z' -J colors of 5.5 < z < 6.7 quasars have a minimum 
atz = 5.8 (Willott et al. 2009). 

4. CFHQS AND SDSS Z = 6 QUASAR SURVEYS 

4.1. CFHQS 

The sample of CFHQS quasars used to derive the luminos- 
ity function in this paper contains all the quasars discovered 
from patches (contiguous sub-regions) which have complete 
near-IR imaging and spectroscopic follow-up to the full mag- 
nitude and color limits. The CFHTLS Wide areas are not in- 
cluded because none of the Wide patches are yet complete. 
All the quasars from the CFHTLS Deep/SXDS, RCS-2 and 
CFHTLS VW are located in complete patches. Therefore the 
CFHQS sample in this paper consists of 12 quasars from the 
RCS-2 (including CFHQS J0050H-3445), 4 quasars from the 
CFHTLS VW and 1 quasai- from the CFHTLS Deep/SXDS. 
The effective sky areas are 494 sq. deg. for the combina- 
tion of RCS-2 and CFHTLS VW and 4.47 sq. deg. for the 
CFHTLS Deep/SXDS. We use the fifth data release (T0005) 
of the CFHTLS Deep. One of the quasars included in the 
CFHQS RCS-2 sample was actually discovered by the SDSS 
(SDSS J23 15-0023 at z = 6.12; Jiang et al. 2008), but as ex- 
plained in Willott et al. (2009) it was already in our follow-up 
list before we learned of its discovery. Since this quasar lies 
within the CFHQS area, it must be included in our sample if 
we are not to underestimate the number of quasars. 
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Figure 3. Left: Photometric completeness as a function of artificial source magnitude and color for a typical RCS-2 patch, CDE2329. White corresponds to 
fractional completeness of 1 and black is 0. The dotted lines show the selection criteria z' < 22.5 and i — z> 2. Note that a small fraction of objects with true 
magnitudes and colors outside of the selection criteria are selected due to photometric errors. Right: Completeness as a function of redshift for the same patch 
using the set of cloned lower redshift quasars. Lines are plotted at three different apparent magnitude levels. 



The patches have variable optical imaging depth and there- 
fore careful completeness corrections have to be applied. The 
method to determine completeness consisted of inserting 1.5 
million artificial point sources into our data and determining 
the fraction recovered by our selection pipeline as a function 
of magnitude and color. Extensive details of this process are 
given in the CFBDS paper of Reyle et al. (2009). An ex- 
ample of the resulting completeness as a function of apparent 
magnitude and color is shown in the left hand panel of Fig. [3] 

For the CFHQS quasars, rather than the brown dwarfs in the 
CFBDS, several further steps are necessary. Galactic extinc- 
tion affects both the observed color and absolute magnitude, 
so is corrected for based on the maps of Schlegel et al. (1998). 
For quasars we need to determine how the color selection 
cut and completeness affects the redshift completeness. This 
is particularly important at z « 5.8 where the typical quasar 
color is approximately the same as our color lower limit. We 
use the 180 cloned z « 3.1 SDSS quasars first discussed in 
Willott et al. (2005a) to determine how the completeness as 
a function of color translates to completeness as a function of 
redshift. This process assumes that the spectra of quasars do 
not evolve significantly from z = 3 to z = 6 in accord with most 
observations (Fan et al. 2004). The only factor that evolves 
with redshift is the transparency of the IGM to Ly a radiation. 
This is modeled using the empirical fit of Songaila (2004). 
The proximity of the CFHQS quasars in color-color space to 
their expected locations (Willott et al. 2009) shows that the 
cloning process is very effective. 

The right hand panel of Fig. [3] shows the completeness as 
a function of redshift for a typical patch from the RCS-2. 
Curves are plotted for three different apparent magnitudes. 
This shows that the completeness does not change much from 
z' = 21' to z' = 22. By the time the magnitude limit is reached, 
the completeness for this RCS-2 patch is about half of its peak 
value. Part of this decline is due to some regions having 10 a 
z' limits brighter than z' = 22.5 and part is due to photomet- 
ric scatter (the hard cut at z' = 22.5 misses some objects with 
intrinsic magnitudes of z' = 22.4). The /'-z' > 2 cut selects 



quasars atz?J5.8. Most of the lowest redshift CFHQS quasars 
are expected to be close to the magnitude limit, because the 
higher magnitude errors increase the number of objects with 
true colors i' -z' < 2 appearing in the selection box. AtzJ$6.5 
quasars move out of the selection box because their z' -J col- 
ors become like those of late T dwarfs and the absolute mag- 
nitude at a fixed z' magnitude decreases rapidly. 

The final step is to convert the CFHQS completeness as 
functions of redshift and apparent magnitude into functions 
of redshift and absolute magnitude. This is accomplished by 
using the 180 cloned lower redshift quasars to mimic the ex- 
pected variance in the relationship between apparent and ab- 
solute magnitude at a fixed redshift. Fig.|4] shows two pan- 
els of the completeness as a function of redshift and absolute 
magnitude for CFHQS. One panel is for all the CFHQS RCS- 
2/VW patches and one is for the combination of the CFHTLS 
Deep and SXDS. The locations of quasars in the complete 
samples used in this paper are also marked. Most quasars 
were found close to the edge of the selection function in ab- 
solute magnitude. The two quasars at z > 6.2 in very dark 
regions of the plot are in regions where the completeness is 
only ^ 5%. The completeness is so low here because only 
parts of some patches have high completeness to the magni- 
tude limit of z' = 22.5. Effectively, this is equivalent to only a 
fraction of the full sky area having been surveyed to the full 
depth. 

4.2. SDSS main sample 

The bright end of the z = 6 quasar luminosity function is 
best sampled by the SDSS main sample described in Fan et al. 
(2000; 2001; 2003; 2004; 2006b). The most recent pubHshed 
sample from Fan et al. (2006b) contains 14 quasars at 5.74 < 
z < 6.42 with uniform selection criteria found from a total 
sky area of ~ 6600 sq. deg. The quasars are selected to a 
magnitude limit of z' < 20.2 and color cut of i' -z' > 2.2. 
Note that due to differences between the filters used by SDSS 
and CFHT, the SDSS sample has better completeness at z ~ 
5.8 than the CFHQS, despite the apparently higher color cut 
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Figure 4. Completeness as a function of absolute magnitude and redshift for 
all SDSS and CFHQS samples used in this paper. White equals peak com- 
pleteness and black equals zero. The red circles show the quasars discovered 
by each sample. SDSS deep z! < 21 corresponds to Jiang et al. (2008) and 
SDSS deep z' > 21 to Jiang et al. (2009). 

of the SDSS. This is the reason why the CFHQS sample has 
a slightly higher median redshift (z = 6.05) than the SDSS 
sample (z = 6.0). 

The SDSS main sample completeness as a function of abso- 
lute magnitude and redshift was supplied by X. Fan and is that 
shown in Fig. 7 of Fan et al. (2003). Fan et al. use a slightly 
different cosmological model to us and therefore their selec- 
tion function and quasar absolute magnitudes are both con- 
verted to our cosmology. 

The SDSS main sample completeness as a function of abso- 
lute magnitude and redshift is also shown in Fig. HI The SDSS 
quasars cluster very strongly to the edge of the selection func- 
tion indicating a fairly steep bright end slope to the luminosity 
function. Indeed, Fan et al. (2004) derived a slope of [3 = -3.2 
from these data. Two SDSS main quasars are found in regions 
of very low selection probability (^ 5%), as was also found 
for the CFHQS. 

4.3. SDSS deep stripe 

A small portion of the SDSS, the deep stripe, has repeated 
imaging which allows one to search for fainter quasars than 
can be found in the main sample. Jiang et al. (2008) con- 
structed a sample of six quasars at 20 < z' < 21 in a sky area 
of 260 sq. deg. In a following paper they extended their search 
to fainter fluxes and found a complete sample of four quasars 
at 21 < z' < 21.8 in 195 sq. deg (Jiang et al. 2009). These 10 
quasars have redshifts between 5.78 and 6.12. 



The completeness for the SDSS deep stripe has been deter- 
mined by Jiang et al. (2008; 2009). In their first paper, quasars 
were selected using follow-up in either the J or H filters (the J 
filter is optimum, but was unavailable for some observations). 
Therefore there are two different selection functions for the J 
and H samples (Fig. 5 of Jiang et al. 2008). We combine the 
J and H selection functions in proportion to the fraction of 
quasars discovered with each filter. Note that there are only 
small differences between these two functions and mostly at 
the high redshift end where no quasars were detected. The 
deeper sample of Jiang et al. (2009) covers some of the same 
area as their brighter sample. We therefore used the brighter 
sample selection function for 65 sq. deg. and for the remain- 
ing 195 sq deg. set the selection function to be the greater of 
either the bright or faint sample functions. This was necessary 
because at some locations in the Mi 450,2 plane, the brighter 
selection function is actually higher than the fainter selection 
function. 

Fig.m shows two panels for the completeness of the two 
SDSS deep stripe samples. The SDSS deep z' > 21 sample 
reaches to within half a magnitude of the CFHQS RCS-2/VW 
sample, but covers less than half the area on the sky. It is 
apparent that quasars from both the SDSS deep samples do 
not cluster as strongly to the selection function as either the 
SDSS main or CFHQS quasars. 

As mentioned above, SDSS J23 15-0023 from Jiang et al. 
(2008) was also discovered in CFHQS imaging and is the only 
quasar to have been identified in the small fraction of overlap 
in sky area between the CFHQS and SDSS deep stripe. This 
quasar will be included in both the CFHQS and SDSS deep 
stripe samples for the purposes of calculating the luminos- 
ity function and the sky areas of the samples will be treated 
as if they are independent. Excluding this duplication, the 
combined CFHQS and SDSS z = 6 quasar sample used in this 
paper consists of 40 quasars at redshifts 5.74 < z < 6.42. 

5. Z=6 QUASAR LUMINOSITY FUNCTION 

5.1. Binned luminosity function 

To help determine the appropriate parametric form of the 
luminosity function we first carry out a simple binning of the 
data. As is well known, the main disadvantage of this method 
is that for small samples with inhomogeneous coverage of pa- 
rameter space, the choice of bins can make large differences 
to the results. In addition, binned luminosity functions are dif- 
ficult to extrapolate beyond well-sampled parts of parameter 
space. 

We use the binned 1 /y„ method of Avni & Bahcafl (1980) 

where Vi is the co-moving volume available for a source j 
in a bin with sizes AMi45() and Az. Due to the small num- 
ber of quasars and small range in redshift, we only use one 
redshift bin. In the following section we will discuss the 
likely amount of evolution across this bin. The available vol- 
ume takes into account the selection functions, p(Mi45o,z), 
described in Sec. |4] and can be calculated as 

^" " II ^^^^450, z) ^ dzdMusQ. 

The volume element dV / dz accounts for the sky area of the 
survey(s) at this location in Mi 450, z space. The luminosity 
function, $(Mi45o), is then calculated as 



«'(Mi450) = 



N 

E 



(AM1450)"'. 
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Figure 5. Binned (open diamonds SDSS main, open triangles SDSS deep 
stripe and filled circles CFHQS) and parametric best fit (thick line) z = 6 
luminosity functions. The thin grey lines show the result of 100 bootstrap 
resamples. The two dashed lines are the power law fits of Jiang et al. (2009). 



The absolute magnitude bins were chosen to avoid human- 
induced bias as much as possible. No attempt was made to 
get approximately equal numbers of objects in each bin, since 
this can bias the results for small samples. 

Unlike in the work of Jiang et al. (2009) which treated their 
z' < 21 and z' > 21 samples separately , we combine the two 
SDSS deep stripe samples of Jiang et al. (2008; 2009) so that 
the SDSS deep quasars are in bins determined by their M^so 
rather than their apparent z' magnitude. 

Fig. |5] shows the binned luminosity function derived from 
the CFHQS and SDSS samples. Data points are plotted at the 
mean absolute magnitude of objects within the bin. Separate 
points are plotted for the CFHQS, SDSS main and SDSS deep 
stripe samples, even though there is some overlap in absolute 
magnitude of the three samples. Note that due to a different 
choice of bins sizes, different SDSS main sample size (Jiang 
et al. included some unpublished quasars) and different cos- 
mological parameters, the SDSS points are not in exactly the 
same locations as in Jiang et al. (2009). 

The dashed lines in Fig. |5] are the two power law fits to the 
binned luminosity function by Jiang et al. (2009), corrected 
for our cosmology. The steeper slope line (/3 = -2.9) is a fit to 
the SDSS main and SDSS deep z' < 21 binned points and the 
flatter slope ((3 = -2.6) is a fit to the SDSS main and all SDSS 
deep binned points in Jiang et al. (2009). These power laws 
from Jiang et al. lie above the SDSS binned points we have 
calculated. The likely reasons for this difference are because 
we used different bins and have a different SDSS main sam- 
ple size. We set the bins to cover absolute magnitude ranges 
without considering the distribution of data points within each 
bin. Jiang et al. (2008, 2009) set the edges of their absolute 
magnitude bins to be equal to the maximum and minimum 
absolute magnitudes of the objects within the bins (L. Jiang, 
priv. comm.). For very few objects in the bin this necessarily 
biases the estimate of the volume available in the bin to low 
values and hence unreasonably high space densities. This is 



particularly apparent for the brightest SDSS deep stripe point 
and faintest SDSS main sample point. A power law fit to our 
SDSS bins would show a somewhat flatter slope and lower 
normalisation than either of those found by Jiang et al. How- 
ever, the discussion above highlights the perils of determining 
the luminosity function from binned data and we therefore 
leave any parametric derivation to the following section. 

There are three CFHQS binned points on Fig.|5] The 
faintest point at M^so « -22 contains the single quasar from 
the Deep/SXDS. The other two points contain all the quasars 
from the RCS-2/VW. The point at M1450 w -24.7 contains 
12 quasars and has by far the smallest statistical error of any 
points on the plot. The two RCS-2/VW points lie below the 
luminosity function defined by the SDSS main and deep stripe 
binned data. However, these data are all consistent within the 
size of the 1 <t error bars, so we do not consider the SDSS and 
CFHQS data to be in conflict. Considering all the data points 
together, there appears to be some evidence for a flattening 
of the luminosity function going from high to low luminosity. 
The main evidence for a break at M1450 « -25 is based on the 
two faintest CFHQS bins. Unfortunately, the small sky area 
probed by the Deep/SXDS and the single quasar found means 
that the space density at very low luminosities is not strongly 
constrained. We will revisit this issue in the following section. 

5.2. Parametric luminosity function 

To overcome the limitations of the binned method, we also 
fitted a parametric luminosity function model to the data. As 
shown by the binned data in Fig.|5] the luminosity function 
is fit well by a power law at the high luminosity end, with 
the possibility of a break at low luminosity. In keeping with 
previous work, we therefore consider the quasar luminosity 
function at z = 6 to follow a double power law function: 

. ^ 10^'-^-^> $(Mr45o) 

'i'\mi45Q,Z) ^Q0.4(a+1)(M,450-M*,5„)_,_ ^q0.4(/3+1)(M,450-M*,5„) ' 

We include an evolution term because the quasar space den- 
sity is well known to decline with redshift from a peak at 
z ~ 2.5 (e.g. Richards et al. 2006). Our sample has only a 
small redshift range so we do not attempt to fit for the evo- 
lution parameter k, but adopt the value of ^ = -0.47 which is 
derived from the evolution from z = 3 to z = 6 for the bright 
end of the luminosity function (Fan et al. 2001). This rate of 
evolution corresponds to a factor of 2 decrease in space den- 
sity across our redshift range from z = 5.8 to z = 6.4. Hence it 
is important to include because redshift and absolute magni- 
tude can be correlated in these color- and magnitude-limited 
samples (e.g. Fig.|4]l. 

This leaves four parameters to be determined:the normal- 
isation $(MJ'45,)), the break magnitude M'^^^q, the bright end 
slope (3 and the faint end slope a. As can be seen in Fig.|5] 
there are very few quasars in our sample at magnitudes well 
below the break in the luminosity function. Therefore we have 
decided to fix the faint end slope to a = -1.5 as indicated by 
AGN surveys at lower redshift (Croom et al. 2004; 2009; 
Hunt et al. 2004). The remaining three parameters are deter- 
mined via a maximum likelihood fit (Marshall et al. 1983). 
The aim is to minimise the function S which equals -21n£, 
where C is the likelihood: 



S = 



N 



1=1 



In [$(Mi45o;,Z,) p{M 1450 i,zd] 
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Figure 6. The contours show the confidence regions (68.3%, 90%, 95. 4%, 
99%) for the bright end slope, /3, and the break absolute magnitude, M|*^jq 
at z = 6. The best-fit z = 6 parameters are marked with a star. The pink, red 
and orange symbols show data at lower redshifts. The pink triangle is z = 3.2 
from Siana et al. (2008) based on SWIRE+SDSS. The orange squares joined 
by a line are from the evolving /3 model fit to the SDSS by Richards et al. 
(2006) from 7, = 2.5 (right) \o z = 4.5 (left). The red dots joined by lines are 
narrow redshift slices from z = 0.75 to z = 2.25 based on the 2SLAQ-I-SDSS 
work of Croom et al. (2009). 

+2 jj -I>(Mi45o,z) p(Musi.),z) ^ dzdMuiQ 

where the first term is a sum over each quasar and the sec- 
ond is integrated over the full possible range of redshift and 
absolute magnitude. Initial starting values of the free parame- 
ters were estimated from the binned luminosity function. The 
parameters were optimised using the amoeba routine (Press et 
al. 1992). 

The best-fit luminosity function parameters are ^(M^^^q) = 
1.14 X 10-'^Mpc"Vag-i, M*45o = -25.13, (i = -2.81. This 
function is plotted on Fig.|5]as the thick line. It passes through 
the 1 a error bars on each of the binned data points showing 
good agreement between the maximum likelihood and binned 
methods. The z = 6 quasar space density at M1450 = -25 is a 
factor of two smaller than that of the /3 = -2.6 fit to binned 
data of Jiang et al. (2009). This has important consequences 
for the quasar ionizing photon output at z = 6 as we shall see 
later. 

Uncertainties on the luminosity function parameters were 
determined assuming a distribution of AS (= S- Smin) 
(Lampton, Margon & Bowyer 1976). Due to the fact that the 
parameters are highly correlated the most useful constraints 
to consider are those on the joint probability of break abso- 
lute magnitude and bright end slope. Note that the faint end 
slope was always fixed at -1.5, so any extra uncertainty due 
to this is unaccounted for 

Fig.|6]plots the x' confidence regions for the combination 
of break absolute magnitude and bright end slope. The pa- 
rameters are obviously correlated such that bright break mag- 
nitudes are consistent with steep bright end slopes. The corre- 
lation is easy to understand because a bright break combined 
with a steep bright end slope would give a similar ratio of 



bright SDSS to faint CFHQS quasars as a luminosity function 
with a faint break and flatter bright end slope. This correla- 
tion tends to disappear when the break magnitude is fainter 
than MJ'45Q w -25 because we have few quasars at such faint 
levels. The 68.3% confidence region extends off the bottom of 
this plot, which effectively corresponds to an unbroken power 
law down to our faintest quasar (M1450 = -22.2). 

The reliability of the x^ parameter uncertainties were 
checked by bootstrap resampling of the data and re-running 
the maximum likelihood fit on these samples (Press et al. 
1992). 100 bootstrap trials were performed. A comparison 
of the fraction of bootstrap trials within various regions of the 
confidence contours in Fig.|6] showed very good agreement. 
25% of the bootstrap trials gave M*^^q > -22 with /3 in the 
range -2.67 to -2.27, which effectively means that the power 
law is unbroken down to our faintest quasar. Although for- 
mally allowed by the data, as mentioned above, this would 
be quite unusual given the double-power law form and break 
magnitudes observed at lower redshifts. The 100 bootstrap 
trials are shown in thin grey lines on Fig.|5] to illustrate the 
range of possible luminosity functions allowed by the data. 
The space density is very well constrained over the range 
-27.5 < M1450 < -24.7 with the 1 a uncertainty being less 
than 0. 1 dex. 

In order to try to break the degeneracy between break mag- 
nitude and bright end slope, and to assess evolution of the 
quasar luminosity function, we consider the results from stud- 
ies at lower redshift. Croom et al. (2009) presented the quasar 
luminosity function at 0.4 < z < 2.6 based on the 2SLAQ and 
SDSS surveys. They had sufficient spread in luminosity and 
redshift to be able to determine the luminosity function in nar- 
row (Az = 0.1) redshift slices. They found evolution in the 
break magnitude and steep end slope (note that a and f3 are 
swapped in Croom et al. compared to our definition). The 
evolution in /3 is quite mild with a steepening towards higher 
redshift. Most of the evolution is in the break magnitude 
which brightens by 4 magnitudes from z = to z = 2 (this is the 
so-called pure luminosity evolution that characterised early 
derivations of the evolving quasar luminosity function, e.g. 
Mathez 1976). The break magnitude and (3 found by Croom et 
al. (2009) are plotted on Fig.|6](we converted from Mg{z = 2) 
to Ml 450 using the relation in Croom et al. and adapted for our 
cosmology). The parameters at 0.75 < z < 2.25 are always 
located at regions in the plot outside of the z = 6 95% confi- 
dence region, with the low redshift slope being steeper than 
that at z = 6 for a given break magnitude. 

Richards et al. (2006) determined the quasar luminosity 
function up to z = 5 using the SDSS DR3. At z > 2.4 they fit 
a model with a redshift-dependent bright end slope and fixed 
break magnitude (due to the fact they did not have any quasars 
below the break at high redshift). Their fit showed the bright 
end slope flattening from (3 = -3.1 at z = 2.4 to /3 = -2.5 at 
z = 4.5. This is shown on Fig.|6]by the orange line with points 
at z = 2.5,3.5,4.5. Support for the choice of break magni- 
tude adopted by Richards et al. (2006) comes from the lu- 
minosity function derived by Siana et al. (2008) using SDSS 
and SWIRE. In fact, the break magnitude and slope found by 
Siana et al. (2008) for z = 3.2 (pink triangle on Fig.|6]l are 
very close to the minimum x^ best fit at z = 6. This would 
suggest little evolution in these parameters from z = 3 to z = 6, 
however, there is a considerable range at 95% confidence for 
z = 6. 

The studies at lower redshift have shown that at all red- 
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shifts in the range 1 < z < 4, the break magnitude lies between 
—26 < Mi45o < -24. Therefore we assume that this also holds 
true at z = 6 in order to provide tighter constraints on /3. Under 
this assumption we find the 95% confidence ranges for (3 to 
be -2.9 < /3 < -2.3 for M*45o = -24 and -3.8 < /? < -2.7 at 

Throughout all the above analysis, the faint end slope was 
fixed at a = -1.5 based on evidence for a at lower redshift 
(Croom et al. 2004; 2009; Hunt et al. 2004). We have checked 
to see how our results would change if instead we adopted 
a different value for a. We re-ran the maximum likelihood 
fitting assuming a = -1.8 and determined best-fit parameters 
of HMU50) = 2.55 X 10-''Mpc-^mag-i, M*^^^ = -26.39 and 
f] = -3.26. Although these best-fit parameters appear sub- 
stantially different to those for a = -1.5, the two luminos- 
ity functions are very similar over the luminosity range cov- 
ered by the SDSS and CFHQS quasars. The steeper bright 
end slope and higher break luminosity are shifted along the 
same direction as the correlation in Fig.|6] The best-fit break 
at M*45o = -26.39 is more luminous than as measured at lower 
redshifts and as noted above we find it unlikely that the break 
would be this luminous at z = 6. Future large area, deep z ^ 6 
quasar surveys are necessary to determine the faint end slope. 

5.3. Consistency of luminosity function with redshift 
distributions 

These colour-selected quasar samples can only find quasars 
within a certain redshift range as shown in Figs.[3]and|4] The 
successful derivation of a luminosity function from these data 
depends upon accurately determining how the completeness 
depends upon redshift and luminosity. In order to check the 
consistency of the observed redshift distributions with the se- 
lection function we have used the selection functions and the 
best-fit luminosity function to generate expected redshift dis- 
tributions for 3 samples: SDSS main, SDSS deep stripe (com- 
bined z' < 21 and z' > 21) and CFHQS RCS-2A'W. These 
expected redshift distributions were then compared to the ob- 
served redshift distributions. A Kolmogorov-Smirnov test 
was applied to determine the probability that the redshift dis- 
tributions were consistent. The results showed that all sam- 
ples are consistent with their expected redshift distributions 
(at probabilities 0.48, 0.78, 0.27, respectively). Jiang et al. 
(2009) noted that their SDSS deep z' > 21 sample had no 
quasars at z > 6.1 and showed that this lack of high-redshift 
quasars was not statistically significant. Our findings agree 
and show that all the samples have redshift distributions as 
expected from the colour selection criteria. 

6. APPLICATIONS OF THE LUMINOSITY FUNCTION 

6. 1. The quasar intergalactic ionizing flux at z = 6 

The constraints placed on the low luminosity end of the 
z = 6 quasar luminosity function here are considerably tighter 
than in previous works (Willott et al. 2005a; Shankar & 
Mathur 2007; Jiang et al. 2009). Therefore it is useful to 
consider the ionizing photon output of the quasar population 
at z = 6. We use the 100 bootstrap trial luminosity func- 
tions to determine the plausible statistical range of the pho- 
ton output. Since the luminosity function is not constrained 
at M*4gQ > -22, we set all the bootstrap models with break 
magnitudes fainter than this to have a break to a = -1.5 at 
M*4g,) = -22. Note however that we do not include ranges on 
other important parameters, such as the faint end slope a (as- 
sumed -1.5 here), the quasar UV spectral slope (assumed that 
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Figure 7. The probability di.stribution for the ionizing photon emission rate 
density from quasars at z = 6. The distribution comes from the range of possi- 
ble luminosity functions consistent with the data determined by the bootstrap 
method. The dashed line shows the required photon emission rate density to 
balance combinations and keep the universe ionized. 

of Telfer et al. 2002) or the photon escape fraction (assumed 
100%). 

The resulting probability distribution in ionizing photon 
rate density is shown in Fig.|2l The peak of the distribu- 
tion is at 4 X 10"*^ photons s"' Mpc""* and the median is at 

5 X 10"*^ photons s"' Mpc"-'. The probability distribution is 
bimodal because the group of bootstrap luminosity func- 
tions with M^^^^1 > -22 have a large number of faint quasars 
and total emission rate density of > 10**^ photons s"' Mpc""*. 
Our results are significantly lower than that previously as- 
sumed of 2 X 10"*^ photons s"' Mpc"^ (Meiksin 2005; Bolton 

6 Haehnelt 2007). Because low-luminosity quasars could 
dominate the ionizing flux, the bootstrap resampling was also 
run for a steeper faint end slope of a = -1.8 and in this case 
the probability distribution is only shifted 0. 1 dex higher than 
for a = -1.5. 

Also, plotted on Fig.|7]is a dashed line to denote the photon 
rate density required in order to balance recombinations and 
keep the universe ionized. We follow Meiksin (2005) to calcu- 
late this number and assume that the dumpiness factor C = 5. 
It is evident from this plot that the quasar population, even 
including the possibility of a large number of faint quasars, is 
insufficient to get even close to the required photon emission 
rate density. Our estimate of the photon rate density is be- 
tween 20 and 100 times lower than the required rate. This is 
consistant with the constraints placed on this rate by the unre- 
solved X-ray background (Djikstra et al. 2004). In compari- 
son, the known galaxy luminosity function at z = 6 shows that 
galaxies provide between 10^" and 4 x 10^" photons s"' Mpc"-' 
(Ouchi et al. 2009). 

6.2. Search strategies for quasars at even higher redshifts 

Now that ^ 50 quasars are known at redshifts 5.7 < z < 6.5, 
the next challenge is to identify quasars at even higher red- 
shifts. This is particularly important given the possible rapid 
evolution in the ionization state of the IGM at these redshifts 
(e.g. Fan et al 2006c). As has been well documented, it is 
challenging to discover quasars at higher redshifts as the Ly a 
line and continuum move into the near-IR. Several large sur- 
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Figure 8. Predictions for z > 6.5 quasar surveys based on our luminosity function. Each plot shows, for a particular redshift slice, the sky surface density of 
quasars brighter than the H magnitude limit. The grey lines show the 100 bootstrapped luminosity functions consistent with our data evolved to higher redshift as 
described in the text. The symbols show the depths of near-IR surveys and the sky density necessary for them to detect one quasar in the redshift slice (UKIDSS 
- red, VISTA - blue, Euclid/JDEM - green). 



veys are now underway or planned aiming to discover quasars 
at z > 7. Given our latest determination of the luminosity 
function at z = 6, we can predict the number of quasars that 
may be detectable at higher redshifts. The main aim of these 
predictions is to investigate the optimum observing strategy 
in a fixed observing time, i.e. wide-shallow vs narrow-deep. 

We calculate the density of quasars on the sky as a func- 
tion of limiting apparent H magnitude in three redshift slices 
(6.5 <z< 7.5, 7.5 < z < 8.5, 8.5 < z < 9.5). The calcu- 
lations are performed for all 100 of the bootstrap resampled 
luminosity functions. The evolution is assumed to continue as 
density evolution at the same exponential rate as determined 
for luminous quasars by Fan et al (2001) from z = 3 to z = 6 
(i.e. 10**-'"^' with k = -0.47). The actual quasar evolution is 
obviously completely unknown and it is possible that some 
luminosity-dependence of the evolution occurs over this red- 
shift range. Theoretical models based on Eddington-limited 
accretion of the z = 6 quasar population traced back to ear- 
lier epochs predict that the black hole accretion rate density is 
somewhere between 10 and 100 times lower at z = 9 than z = 6 
(Li et al. 2007; Sijacki et al. 2009). Note that these models 
are focussed on overdense regions, due to simulation limita- 
tions, so may not represent typical regions. However, their 
evolution compares well with our factor of ^ 30 evolution in 
space density from z = 9 to z = 6. 

Recent deep imaging from WFC3 onboard Hubble Space 
Telescope and from the ground has enabled the identification 
of galaxies up to z = 8.5 (Ouchi et al. 2009; Bouwens et al. 
2009; Oesch et al. 2009). Assuming these galaxies are not 
mostly low redshift interlopers (see Capak et al. 2009), the 
evolution in the galaxy luminosity function from z = 6 to z = 7 
and z = 8 is not too rapid. Bouwens et al. (2009) and McLure 
et al. (2009) show the faint end evolution from z = 6 to z = 8 
is less than a factor of ten. Furthermore, Labbe et al. (2009) 
show the stellar mass density evolves only by a factor of three 
between z = 6 and z = 7. Although it is unclear how the evo- 
lution of mostly low luminosity galaxies relates to that of lu- 
minous quasars, these results are consistent with the evolution 
for quasars we have adopted. 

Fig. [8] shows the 100 bootstrap quasar count predictions 
in each of the three redshift slices. The expected counts 
are most strongly constrained for quasars with similar ab- 



solute magnitude to the SDSS and CFHQS quasars, corre- 
sponding to observed magnitudes 20 to 23. Fontanot et al. 
(2007a) made similar predictions for the near-IR quasar num- 
ber counts based on extrapolation of the luminosity functions 
of Fontanot et al. (2007b) at z 4 and Shankar & Mathur 
(2007) at z ^ 6. They did not show results for the bright end 
of the luminosity function at apparent magnitudes H < 24. At 
H > 24- our results are consistent with the predictions based 
on the evolving model of Fontanot et al. (2007b). 

We now consider the specifics of various ongoing and fu- 
ture near-IR sky surveys which might be able to discover high- 
z quasars. All limiting magnitudes are given as IOct AB mag- 
nitudes, because experience of the CFHQS shows that this 
level of photometric accuracy in the bands redward of the 
Ly a line is necessary for reliable colour selection. It should 
also be noted that we simply consider whether these surveys 
cover sufficient depth and area at H band to be able to de- 
tect the quasars. Successful quasar discovery depends upon 
colour criteria that isolate the quasars from all contaminants 
such as brown dwarfs. This usually means that data in filters 
shortward of the Ly a break must go 1 to 2 magnitudes deeper 
than the H band data. 

The UKIRT Infrared Deep Sky Survey (UKIDSS; 
Lawrence et al. 2007) began in 2005 and should be com- 
pleted in a few years. It has several components of which the 
two we consider here are the Large Area Survey (LAS; 4000 
sq. deg. toH= 19.5) and the Ultra Deep Survey (UDS; 0.77 
sq. deg. to H = 24.5). A primary goal of the LAS is to find 
z = 7 quasars (Hewett et al. 2006) and it so far has found 
two at z « 6 (Venemans et al. 2007; Mortlock et al. 2009). 
In Fig. [8] these surveys are plotted alongside the 100 bootstrap 
predictions such that the location on the plot indicates that one 
quasar would be found in the survey in that redshift slice. At 
z 7 the LAS falls below the predictions indicating that be- 
tween 1 and 3 quasars would be expected in this survey at this 
redshift. The UDS lies right at the upper edge of the predicted 
quasar counts meaning that it is likely that it will contain < 1 
quasar at 6.5 < z < 7.5. Note that the UDS forms part of the 
deep component of the CFHQS and the lone deep CFHQS 
quasar at z = 6.01 actually lies within the UDS. At higher red- 
shifts, the LAS may contain ~ 1 quasar at 7.5 < z < 8.5, but 
it is very unlikely for the UDS to contain any higher redsift 
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quasars. In conclusion, the wide area component of UKIDSS 
is much better suited to discovering high-z quasars than the 
deep component. 

Set to commence imminently, the ESO VISTA Telescope 
surveys (Sutherland 2009) are more efficient than UKIDSS 
due to the use of a larger camera. There are four compo- 
nents of the VISTA surveys that could potentially detect high- 
z quasars. These are VISTA Hemisphere Survey (VHS; 5000 
sq. deg. to H = 19.9 coincident with the Dark Energy Survey 
in the optical), VISTA Kilo-Degree Infrared Galaxy Survey 
(VIKING; 1500 sq. deg. to H = 19.9 coincident with VST 
KIDS in the optical), VISTA Deep Extragalactic Observations 
Survey (VIDEO; 15 sq. deg. to H = 23.7) and UltraVISTA 
(0.73 sq. deg. to H = 25.4 in the COSMOS field). As shown 
on Fig.|8]these surveys should be deep enough at H to detect 
many quasars at 6.5 < z < 7.5 and a few at 7.5 < z < 8.5. 
The most suitable surveys are the wider area surveys of VHS 
and VIKING, rather than VIDEO and UltraVISTA. This is a 
consequence of the relatively shallow z = 6 luminosity func- 
tion we have found. With a slope of (3 > -3.0, more quasars 
can be found in a wide area shallow survey than a narrow area 
deep survey in a fixed observing time. However, we caution 
that if the bright end of the luminosity function steepens at 
high redshift, due for example to the lack of time available for 
the build-up of the required mass black holes, then this would 
reverse the situation. Therefore VIKING is probably more 
likely to yield z ^ ^ quasars than VHS. Finally, we note that 
none of the VISTA surveys are likely to detect z ~ 9 quasars. 

What about the longer-term future? Several ambitious dark 
energy projects are being planned such as Euclid and JDEM. 
These would use some combination of supernovae, baryonic 
acoustics oscillations and weak lensing to determine the pa- 
rameters governing the evolution of dark energy. In order to 
make these measurements it is necessary to make deep sur- 
veys of large sky areas in the near-IR and therefore such sur- 
veys would also be useful to detect the most distant quasars. 
Given the early planning stage it is not known yet how much 
sky area will be covered to what depth. For the purposes of 
this exercise we use a preliminary Euclid plan which include 
a wide 20000 sq. deg. imaging survey to H = 23.3 for weak 
lensing and a deep imaging survey of 30 sq. deg. to H = 25.3 
for supernovae. The green circles on Fig.[8]show that the wide 
survey is much more effective at detecting distant quasars due 
to the flatter faint end slope (note this is not strictly a fair com- 
parison because to these depths and areas in a single filter, the 
wide survey requires 17 times as long as the deep survey). The 
wide survey would be able to detect ^ 1000 quasars at z '--^ 7, 

400 at z '--^ 8 and ^ 100 at z ~ 9. As before, we caution that 
sufficiently deep data at filters below the Ly a break are nec- 
essary to provide accurate colours, so this is an upper limit to 
the number of quasars that might be discoverable. We do not 
consider even deeper surveys over small sky areas with JWST, 
because the z = 6 luminosity function is not well constrained 
for such low luminosity quasars. 

7. CONCLUSIONS 

We have presented discovery data for nine new quasars in 
the CFHQS, bringing the total number so far to 19. The 
CFHQS is nearing completion and ongoing follow-up may 
find a few more CFHQS quasars. Further observations of 
CFHQS quasars are being pursued to understand reionization, 
black hole mass growth and host galaxy evolution. 

The CFHQS and SDSS surveys were combined to derive 
the z = 6 quasar luminosity function from a sample of 40 



quasars. The normalisation of the luminosity function is 
found to be « 40% lower than that determined previously us- 
ing binned data (Jiang et al. 2008; 2009). The luminosity 
function is well constrained down to M1450 « -24 but has a 
large uncertainty at M1450 « -22 where our faintest quasar is. 
Future larger sky area, deep surveys, such as VISTA VIDEO, 
are necessary to constrain this part of the z = 6 quasar lumi- 
nosity function. We are also working on measuring black hole 
masses for z > 6 quasars. In a future paper we will combine 
our knowledge of the luminosity function and black hole mass 
distribution to consider the growth of black holes at this early 
epoch. 

We have calculated the quasar ionizing flux based on pos- 
sible realizations of the z = 6 quasar luminosity function and 
found that the ionizing photon output of quasars is between 
20 and 100 times lower than the rate required to keep the uni- 
verse ionized. Predictions have been made for the number of 
z > 6.5 quasars that might be found by current and future sur- 
veys. Whilst UKIDSS may find a few quasars at z '--^ 7, the 
VISTA surveys should find many more and even some up to 
z 8. Obtaining large samples of z ~ 8 quasars and push- 
ing to even higher redshifts will require even more ambitious 
projects such as Euclid or JDEM. 



Based on observations obtained with 
MegaPrime/MegaCam, a joint project of CFHT and 
CEA/DAPNIA, at the Canada-France-Hawaii Telescope 
(CFHT) which is operated by the National Research Council 
(NRC) of Canada, the Institut National des Sciences de 
rUnivers of the Centre National de la Recherche Scientifique 
(CNRS) of France, and the University of Hawaii. This work 
is based in part on data products produced at TERAPIX 
and the Canadian Astronomy Data Centre as part of the 
Canada-France-Hawaii Telescope Legacy Survey, a collab- 
orative project of NRC and CNRS. Based on observations 
obtained at the Gemini Observatory, which is operated by 
the Association of Universities for Research in Astronomy, 
Inc., under a cooperative agreement with the NSF on behalf 
of the Gemini partnership: the National Science Foundation 
(United States), the Particle Physics and Astronomy Research 
Council (United Kingdom), the National Research Council 
(Canada), CONICYT (Chile), the AustraHan Research Coun- 
cil (Australia), CNPq (Brazil) and CONICET (Argentina). 
This paper uses data from Gemini programs GS-2008B-Q- 
23, GS-2009A-Q-3, GS-2009B-Q-25, GN-2008B-Q-42 and 
GN-2009A-Q-2. Some of the data presented herein were 
obtained at the W.M. Keck Observatory, which is operated 
as a scientific partnership among the California Institute of 
Technology, the University of California and the National 
Aeronautics and Space Administration. The Observatory 
was made possible by the generous financial support of 
the W.M. Keck Foundation. Based on observations made 
with the ESO New Technology Telescope at the La Silla 
Observatory. Based on observations made with the Nordic 
Optical Telescope, operated on the island of La Palma jointly 
by Denmark, Finland, Iceland, Norway, and Sweden, in the 
Spanish Observatorio del Roque de los Muchachos of the 
Instituto de Astrofisica de Canarias. Thanks to Howard Yee 
and the rest of the RCS2 team for sharing their data and to 
the queue observers at CFHT and Gemini who obtained data 
for this project. Thanks to Xiaohui Fan, Linhua Jiang and 
Scott Croom for providing information relating to their work 
on the luminosity function and Jasper Wall for interesting 
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APPENDIX 
FINDING CHARTS 

Fig. |9] presents 3' x 3' finding charts for the CFHQS quasars. All images are centered on the quasars and have the same 
orientation on the sky. These are the CFHT MegaCam z'-band images in which the quasars were first identified. MegaCam has 
gaps between the CCDs and these data were not dithered so the gaps remain and are evident in 5 images where the quasars lie 
close to the edge of a CCD. For CFHQS J0050+3445, CFHQS J0136+0226, CFHQS J0316-1340 and CFHQS J1059-0906, the 
Megacam data are a single exposure leading to many cosmic rays in the final images. 




CFHQS J1429+5447 
Figure 9. z'-band finding charts for the CFHQS quasars. 



CFHQS J0221-0802 
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Figure 9. (cont.) 



